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ABSTRACT 

The  problem  of  applying  horizontal  and/or  vertical  parity  checks  for 
error  detection  to  actual  burst  channels  is  considered.  It  is  demonstrated 
that  a  single-dimension  parity  check  will  achieve  two  to  four  orders-of- 
magnitude  improvement,  and,  within  the  limit  of  the  data  sample,  both 
dimensions  together  detect  all  errors. 

It  is  proved  that  the  parity  check  at  the  intersection  of  the  horizontal 
and  vertical  checks  can  be  calculated  from  either  set  of  checks  and  that 
errors  can  be  detected  in  the  information  independent  of  those  in  the  non¬ 
intersecting  checks. 
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PERFORMANCE  OF  TWO-DIMENSIONAL  ERROR 
DETECTION  ON  DIGITAL  HF  AND 
TROPOSCATTER  CHANNELS 


INTRODUCTION 

In  previous  work  Kuhn^  analyzed  the  error-detection  performance 
of  two-dimensional  (horizontal/vertical)  binary  parity  check  codes  in  a  random- 
error  channel  and  predicted  the  undetected  block  error  probabilities.  Owing 
to  a  lack  of  burst  statistics,  Kuhn  gave  only  a  qualitative  analysis  of  perform¬ 
ance  in  real  burst  channels.  In  this  paper  the  error-detection  performance 

[2] 

of  this  scheme  shall  be  described  for  the  previously  reported  high-frequency 

[3] 

and  troposcatter 1  channels  in  which  errors  occur  in  bursts. 

STRUCTURE  OF  THE  CODE 

A  two-dimensional  binary  parity  check  code  is  structured  by  taking  a 
sequence  of  m  binary  digits  from  the  information  source  and  calculating  a 
parity  bit.  This  generates  an  (m  +  l)-bit  word.  The  process  is  repeated 
n  times  and  to  the  n  words  is  appended  a  parity  word.  The  parity  word  is 
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calculated  by  using  the  first  bit  of  each  information  word  to  find  the  first  bit 
of  the  parity  word  and  repeating  the  calculation  for  the  m  information  bits/ 
word.  As  proved  in  the  Appendix,  the  last  bit  of  the  parity  word  can  be  cal¬ 
culated  from  the  word  parity  bits  or  from  the  bits  of  the  parity  word.  Be¬ 
cause  of  the  relationship  of  the  bits  in  the  parity  word  to  the  information  bits, 
the  parity  word  is  frequently  called  vertical  or  column  parity.  Similarly, 
the  word  parity  bits  are  frequently  called  horizontal  or  row  parity.  The 
parity  bits  are  calculated  under  an  even-calculation  rule.  This  rule  requires 
that  the  sum  of  all  bits  in  a  row  (column)  be  even,  where  addition  is  modulo 
2  with  no  carry  and  zero  is  an  even  number.  Thus,  in  encoding,  if  the  sum 
of  the  information  bits  in  a  row  (column)  is  odd  the  parity  bit  is  a  "one.  "  In 
decoding,  each  row  (column)  is  checked  to  see  if  the  sum  is  still  even  (parity 
is  satisfied).  If  it  is  not,  an  error  is  declared  as  detected.  It  is  possible, 
however,  that  an  even  number  of  errors  will  occur  in  a  row  (column).  In 
this  case  the  sum  is  not  changed  and  the  errors  are  not  detected  in  the  row 
(column).  Further,  if  odd  parity  were  originally  used  in  the  encoding  and 
decoding  rules  the  performance  would  be  the  same.  The  implementation  of 
this  code  requires  m  +  1  flip-flops.  Of  these,  m  are  used  for  the  informa¬ 
tion  column  calculations,  and  since  rows  are  processed  sequentially,  only  one 
is  needed  for  the  row  calculations  including  the  calculation  on  the  parity  word. 

REVIEW  OF  CHANNEL  DATA 

In  October  of  1965,  field  tests  were  conducted  on  NRD  high-frequency 

circuits  between  Antigua  and  Ascension  Islands  to  record  HF  digital  error 
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patterns  .  Tests  were  conducted  for  10-minute  periods  (runs)  at  2400 
bits/sec  using  Kineplex  TE-216  16- tone  PSK  modems.  A  52-bit  test  message 
was  transmitted  from  Antigua  to  Ascension  and  back  to  Antigua  on  a  looped 
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basis,  where  the  received  message  was  compared  with  a  suitably  delayed 

replica  of  the  test  message.  Bits  not  in  agreement  were  declared  in  error 

[3] 

and  recorded  for  computer  processing.  The  troposcatter  data1'  were  col¬ 
lected  using  Sebit  24B  vestigial  sideband  AM  modems  at  2400  bits/sec  on  a 
link  dominated  by  a  583-mile  troposcatter  hop.  Test  runs  were  90  minutes 
long.  The  data  are  summarized  in  Table  I. 


Table  I 

Size  of  Data  Sample 

Channel 

High- Frequency 
Troposcatter 


Total  Bits 


Bit  Error  Rate  = 


Total  Errors 
Total  Bits 


4.89  x  10 


5.47  x  10 


.8 


5.26  x  10 


3.95  x  10 
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Analysis  of  the  HF  data1  J  indicates  that  the  errors  occur  in  bursts  of 

approximately  5%  error  density  and  3000-bit  average  length.  Troposcatter 

[3] 

bursts1  J  occur  in  densities  ranging  from  10%  to  a  few  of  75%  and  average 
less  than  1000  bits  in  length. 


PERFORMANCE  RESULTS 

For  purposes  of  analysis  it  is  assumed  that  an  all-zero  message  is 
transmitted.  This  message  satisfies  even  parity.  If  the  error  patterns  are 
added  to  the  message  sequence  the  result  is  the  error-pattern  data.  Thus 
the  code  performance  is  evaluated  by  checking  the  error  patterns  for  even 
parity.  The  errors  were  divided  into  blocks  of  (n  +  1)  rows  and  (m  +  1) 
columns  containing  n  rows  and  m  columns  of  information.  Within  each 
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block  having  errors,  rows  and  columns  were  checked  for  parity.  If  any 
parity  check  was  not  satisfied,  a  block  in  error  was  declared  as  detected. 
Since  the  number  of  blocks  with  errors  was  also  known,  the  undetected  error 
rate  could  be  obtained.  Results  were  also  obtained  using  row  or  column 
parity  check  detection  only. 

As  an  example,  the  results  for  m  =  8  and  n  =  7  on  the  tropospheric 
data  are  as  shown  in  Table  n. 


Table  II 

Performance  of  Ilorizontal/Vertical  Parity  Check 


Block  Size 

Total  Blocks 

Blocks  in  Error 

Case  A:  Blocks  Detectable 
Using  Row  Parity  Only 

Case  B:  Blocks  Detectable 
Using  Column  Parity  Only 

Case  C:  Blocks  Detectable 

Using  Two-Dimensional  Parity 

Channel  Block  Error  Rate 

Undetected  Block  Error  Rate: 

Case  A 

Case  B 

Case  C 


9  columns  x  8  rows  =  72  bits 
7,309,865 
27,181 

25,332 

26,984 

27,177 

.0037 

.00025 

.000026 

.00000055 


From  the  table  the  burst  nature  of  the  channel  is  evident  in  that  only 
27,181  of  7.3  million  blocks  had  errors.  Since  data  are  transmitted 
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sequentially  in  words,  the  bursts  occur  across  rows,  leading  to  poorer  row- 

undetected  error  rates.  In  Figures  1  and  2  the  results  are  presented  for 

both  data  channels  for  the  cases  of  transmission  of  8-bit  (m  =  8)  computer 

bytes  and  square  blocks,  respectively.  For  the  HF  data  there  is  always  a 

block  length  beyond  which  there  are  no  undetected  block  errors,  and  for 

m  =  8  most  rows  have  but  one  error,  thus  favoring  row  detection  over 

column  detection.  In  transmission  of  computer  bytes  through  the  tropo- 

scatter  channel  there  are  four  patterns  of  errors  which  sometimes  go 

undetected,  but  square  blocks  are  always  detectable  on  a  block  basis.  Had 

[4] 

the  messages  been  interleaved  prior  to  transmission,  the  square  block 
row  and  column  performance  results  would  have  been  interchanged. 

CONCLUSIONS 

Within  the  size  of  the  data  sample  used  (6  hrs  HF  and  61  hrs  troposcatter) 
the  limits  of  performance  have  not  yet  been  reached.  One  can  only  speculate 
on  how  much  data  need  be  collected  before  a  valid  measure  of  undetected  error 
probability  is  attained,  but  in  any  event  it  appears  that  a  test  program  to 
collect  such  data  would  be  unpractically  long  and  expensive.  The  results  can¬ 
not  be  compared  to  those  of  Kuhn^  ,  since  to  validate  undetected  error  rates 
-14 

of  10  would  require  60,000,000  hrs  of  data  at  2400  bits/sec.  It  can  be  con¬ 
cluded,  however,  that  simple  parity  gives  excellent  error-detection  perform¬ 
ance  at  a  reasonable  cost.  Furthermore,  two  dimensions  of  error  detection 
are  far  superior  to  either  horizontal  or  vertical  alone.  Additionally,  the  use 

of  this  simple  scheme  gives  lower  undetected  error  rates  than  those  pre- 

[5] 

viously  reported1  for  BCH  codes. 
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Figure  1 .  Performance  of  Two-Dimensional  Parity  Check 
Error  Detection  (Computer  Bytes) 
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NUMBER  OF  WORDS  (W)  IN  W  X  W  BLOCK 


Figure  2.  Performance  of  Two-Dimensional  Parity  Check 
Error  Detection  (Square  Block) 
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APPENDIX 


In  this  appendix  a  theorem  which  can  be  used  in  encoding  a  two-dimen¬ 
sional  parity  check  matrix  is  proved.  It  is  demonstrated  that  the  encoding 
rule  must  be  modified  if  the  dimensions  of  the  matrix  are  not  the  same  in 
kind  (odd  or  even)  and  odd  parity  is  used.  A  corollary  to  this  theorem  per¬ 
mits  the  detection  of  errors  in  transmitted  information  independent  of  the 
errors  in  the  one-dimensional  parity  bits. 


Theorem: 


In  a  two-dimensional  parity  check  matrix,  with  n  information  words 
and  m  information  bits/word,  where  each  word  has  a  parity  bit  appended  to 
it  and  a  parity  word  is  appended  to  the  matrix;  the  bit  of  the  parity  word  which 
corresponds  to  the  word  parity  bits  can  be  calculated  from  either  set  of  parity 
bits. 

Proof  (Case  1  -  Even  Parity  Rule): 

Given  a  matrix  of  information  bits  A  with  elements  a.  .  (a  =  0,1), 

i.J  i,j 

(i  =  1,2 . n)  ,  (j  =  1,2 . m).  With  an  even  parity  constraint,  the  word 

parity  bit  a^  is  given  by 


m 


a. 


i,  m+1 


[modulo  2]  ,  i  =  1,  2, . .  .  ,  n  .  (1) 
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The  block  parity  bit  when  calculated  from  the  word  parity  bits  is  given  by 


a 


n+1,  m+1 


n 


I 


i=l 


a. 

i, 


m+1 


[modulo  2] 


n  m 


[modulo  2]  . 


If  the  sum  is  reordered, 


m 


a 


n+1,  m+1 


[modulo  2] 


where  a  „  .  is  given  by 
n+1, 3  &  J 


n 


a 


n+1,  j 


[modulo  2]  ,  j  =  1,  2, . . 


Comparing  Equations  (2)  and  (4)  it  is  found  that  while  the  bit  £ 
was  originally  calculated  from  the  word  parity  bits  the  result  is  the 
if  it  had  been  found  from  the  bits  of  the  parity  w'ord. 


Proof  (Case  2  -  Odd  Parity  Rule): 

Sub-case  2a  (m  and  n  are  even). 

In  sub-case  2a,  Equation  (1)  is  replaced  by 

m 


a. 

1, 


m+1 


1  + 


a.  . 

i.J 


[modulo  2] 


3=1 


1,2... 


(2) 


(3) 


(4) 


,  m  .  (5) 


n+1,  m+1 
same  as 


•  n  .  (6) 
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It  follows  that 


a  i  =  1  + 
n+1,  m+1 


I  \ 


m+1 


i=l 


[modulo  2j  (7) 


i+ if  +  ia- 

i=l  \  j=l 


i,J 


[modulo  2]  (8) 


n  n  m 

1  +  ^  (1)  +  \  ai  j  [modulo  2] 


i=l 


i=l  j=l 


(9) 


n  m 

i  ♦  y  y  a. 


i.j 


i“l  j=l 


m 

1  +  Y  a 


n+1,  j 


[modulo  2]  (10) 


[modulo  2]  (11) 


j=l 


where 


n 

a  =  1  +  )  a.  .  [modulo  2]  ,  j  =  1,2,  ...m.  (12) 

n+1,  j  L:  l,  J 

i=l 


Sub-case  2b  (m  and  n  are  odd). 


In  this  sub-case,  the  derivation  follows  that  of  sub-case  2a  with  the 

exception  that  the  value  a  „  „  equals  one  plus  that  value  of  sub-case  2a. 

v  n+1,  m+1 
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Sub-case  2c  (m  and  n  are  different  in  kind). 


In  sub-case  2c,  calculation  of  a  ,  yields  two  answers  differing 

’  n+1,  m+1  h 

by  one  depending  on  the  bits  used  in  the  calculation.  In  this  sub-case,  the 
calculation  using  an  even  number  of  terms  must  be  increased  by  one. 

This  proof  can  be  extended  to  an  r-dimensional  parity  check  matrix, 
but  such  matrices  are  not  used  in  practice. 

The  above  theorem  leads  to  the  following  corollary. 

Corollary: 

Parity  bit  a  ,  ,  detects  an  odd  number  of  errors  in  the  informa- 

J  n+1,  m+1 

tion  portion  of  the  block  (including  bit  an4^  m+i^  independent  of  errors  in 
the  other  parity  bits. 

This  is  a  direct  result  of  Equations  (3)  and  (10). 
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